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SUMMARY 

The effects of system parameters on the interface condensation rate in a 
laminar jet-induced mixing tank are numerically studied. The physical system 
consists of a partially-filled cylindrical tank with a si ightly-subcooled jet 
discharged from the center of the tank bottom toward the liquid-vapor interface 
which Is at a saturation temperature corresponding to the constant tank pres- 
sure. Liquid is also withdrawn from the outer part of the tank bottom to main- 
tain the constant liquid level. The jet velocity is selected to be low enough 
such that the free surface Is approximately flat. The effect of vapor super- 
heat is assumed to be negligible. Therefore, the Interface condensation rate 
can be determined from the resulting temperature field In the liquid region 
alone. The nondimensional form of the steady-state conservation equations are 
solved by a finite-difference method for various system parameters Including 
liquid height to tank diameter ratio, tank to jet diameter ratio, liquid Inflow 
to outflow area ratio, and a heat leak parameter which characterizes the uni- 
form wall heat-flux. Detailed analyses based on the numerical solutions are 
performed and simplified equations are suggested for the prediction of inter- 
face condensation rate. 


NOMENCLATURE 

A surface area of central jet, outflow, or Interface 
B tank to jet diameter ratio, D/d 

C p specific heat at constant pressure 

D tank diameter 

d jet diameter 

g gravitational acceleration 


‘Member, AIAA. 


H liquid height from the tank bottom 

h c condensation heat-transfer coefficient 

h^ latent heat of condensation 

Ja Jakob number, - T^)/h^g 

k thermal conductivity 

m condensation mass flux 

c 

Nh wall heat-flux parameter, Nh = q w D/k(T s - T\ ) 

p pressure 

Pg equilibrium hydrostatic pressure 

2 

p* dimensionless pressure, (p - p )/pu. 

W J 

0 jet volume flow rate 

q tank wal 1 heat flux 

r radial coordinate measured from the centerline 

r* dimensionless radial coordinate, r/D 

Re. jet Reynolds number, pu.d/y 

J J 

St c condensation Stanton number, h c /pUjCp 

T temperature 

T* dimensionless temperature, (T - T.)/(T - T.) 

J J J 

u axial velocity 

u c condensation-induced velocity (positive value) 

u* dimensionless axial velocity, u/u^ 

u* dimensionless condensation-induced velocity, u /u i 

v radial velocity 

v* dimensionless radial velocity, v/u^ 

v* dimensionless radial velocity at interface 

x axial coordinate measured from the tank bottom 

x* dimensionless axial coordinate, x/D 


2 


Xp potential core length 

n dynamic viscosity 

p liquid density 


average value over the Interface 

Subscripts : 

j evaluated at jet inlet 

out evaluated at outflow location 

s evaluated at liquid-vapor interface 

c evaluated at condensation condition 

v evaluated at vapor condition 


INTRODUCTION 

The pressure control of cryogenic storage tanks in space is one of the 
major technologies being developed in the NASA cryogenic fluids management 
program- The heat transferred through insulated tank walls can result in 
thermal stratification of the cryogenic propellant which increases the tank 
self-pressurization rate. A preferred method of controlling the pressure of 
a cryogenic storage system Is the use of axial jet-induced mixing (refs. 1 
and 2). A slightly subcooled jet can provide mixing of tank liquid and trans- 
fers heat from the Interface into the bulk liquid. Vapor condensation is 
induced and the pressure of the tank is then reduced. Therefore, the interface 
condensation plays the key role in controlling the pressure of a jet-induced 
mixing tank. 

Several studies of steady-state fluid mixing and condensation rate have 
been conducted. Thomas (ref. 3) measured the condensation rate of steam on 
water surfaces mixed by a submerged turbulent jet and found that the conden- 
sation rate was roughly proportional to the jet Reynolds number. Dominick 
(ref. 4) investigated the effects of jet injection angle and jet flow rate on 
the condensation rate in a Freon 113 tank. The interface heat transfer was 
observed to increase as jet injection angle became more normal to the inter- 
face. However, the average heat transfer coefficient at the interface was 
found to be increasing with jet Reynolds number only to a power of 0.73. Sonin 
et al . (ref. 5) measured the steady-state condensation rate of steam in a water 
tank with condensation process dominated by the liquid-side turbulence near 
the Interface. A correlation between the condensation rate and characteristic 
interface turbulent velocity was developed. Hasan and Lin (ref. 6) used a 
finite-difference method to solve time-averaged conservation equations along 
with a k-e turbulence model for the prediction of turbulent velocity. The 
numerical prediction was in good agreement with Sonin's data except for the 
region close to the interface where the diffusion process was dominant. The 
failure of the numerical solution for the near-interface region, as stated in 
reference 6, was probably due to an inappropriate turbulence model or the 
inappropriate turbulence boundary conditions applied at the interface. The 
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effects of flow parameters, such as jet Reynolds number and Prandtl number on 
the Interface condensation rate in a laminar jet-induced mixing tank was 
studied numerically by Lin (ref. 7). The numerical results showed that the 
average condensation heat-transfer coefficient Increased linearly with the jet 
Reynolds number and varied with the liquid Prandtl number to a power of about 
0.97 for 0.85 < Pr < 2.65. 

Besides the flow parameters, the vapor condensation rate in a mixing tank 
is also dependent on the system parameters which affect directly the design of 
the mixing device. Lacking an adequate design data base, any overdesign of the 
mixing device will increase its weight and power consumption. This will lower 
the performance efficiency of the mixing device in reducing tank pressure. 

This paper presents results of a numerical investigation on the effects of sys- 
tem parameters on the interface condensation rate and provides useful informa- 
tion for the design of a mixing device. 

Reference 7 showed that for a small jet to tank diameter ratio, most of 
the vapor condensation occurs in the central part of the interface so that the 
normal -grav i ty solution (flat interface) may be a good approximation to low- 
gravity condition (curved interface). Only laminar flow is considered in the 
present investigation. However, it is expected that the effects of system 
parameters on the interface condensation process will be qualitatively the same 
for both the laminar and the turbulent jets. 

Numerical solutions of the nondimensional continuity, momentum and energy 
equations are obtained by a finite-difference method. Calculations cover a 
wide range of system parameters such as liquid height to tank diameter ratio, 
tank to jet diameter ratio, liquid inflow to outflow area ratio, and a param- 
eter Nh which characterizes the uniform wall heat-flux. Detailed analyses 
are performed and simplified equations are suggested to predict the effects of 
these parameters on the condensation Stanton number which characterizes the 
steady-state interface condensation rate. 


PROBLEM FORMULATION AND DESCRIPTION 

The physical system and the coordinates used to analyze the problem are 
shown in figure 1. A laminar axial-jet of diameter d is located at the cen- 
ter of the bottom of a cylindrical tank of diameter D. Liquid is continuously 
withdrawn from the outer portion of the tank bottom to maintain a constant liq- 
uid fill level. The tank outflow to jet inflow area ratio (A ou t/Aj) is varied 
from 7.5 to 360 and the tank to jet diameter ratio (D/d) is varied from 12.5 to 
33.3. Both the central jet and outflow velocities, Uj and u ou t, respec- 
tively, are assumed uniform. The effect of vapor superheat is neglected which 
implies that Cp V (T v - T s )hfg << 1. Also, the effect of vapor motion on the 
condensation process is assumed to be negligible. The above assumptions effec- 
tively decouple the liquid region from the vapor region. Therefore, the con- 
densation rate can be determined from the solution of the temperature field in 
the liquid region only. 

Various fill levels are considered with the liquid height to tank diameter 
ratio (H/D) ranging from 0.25 to 1.25. The liquid-vapor interface is assumed 
to be flat (wave free), shear-free and at a constant saturation temperature, 

T s . The axial jet is at a constant subcooled temperature, Tj, and the tank 
outflow region is assumed to have zero temperature gradient. Tank walls are 
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subjected to a uniform heat leak. The dimensionless wall heat-flux parameter, 
Nh (Nh E q w D/k(T s - TO), which characterizes the uniform wall heat-flux is 
varied from 0 to 4. All the associated thermodynamics and transport properties 
are assumed constant since only small temperature variation is considered in 
the whole flowfield. 


The energy and mass balance in the liquid region can be expressed, respec- 
tively, as 


I pU out ( 'p T out c * A out 

pu j c p T j dA j ■ 

J m c h fg dA s * J pu c c pV\ * DHd w 

(1) 


pu out A out - J 

m c dA s . pUjAj 

(2) 


where m c , hf g , C p , T ou t, and q w are the local interfacial condensation mass 
flux, latent heat of condensation, specific heat at constant pressure, local 
outflow temperature, and wall heat flux, respectively. The condensation- 
induced velocity is obtained by 


u c = r 


(3) 


By using equations (2) and (3), equation (1) becomes 


m c dA s 


! J fr^) - " H 


1 + Ja 1 - 


T out - V 

T s - T i . 


(4) 


where the jet volume flow rate (Q) , the Jakob number (Ja), and the average 
outflow temperature are defined as Q = u j A j , Ja = C p (T s - Tj)/hfg, and 
Tout = S T ou tdA ou t/A ou t, respectively. The form of equation (4) is convenient 
to estimate the interfacial condensation mass flux by using the measurable 
experimental variables. The condensation mass flux at the Interface is given 
by 


■cVi 


■ J ”c h fg dA s 


V s 


(5) 


where m c is the average condensation mass flux at the interface, and k is 
the thermal conductivity. Thus, the condensation-induced velocity can be rela- 
ted to the interfacial temperature gradient by 
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Equations (4) and (5) can be used as the termination criterion for the itera- 
tive calculation or as a check of the numerical solutions. 


MATHEMATICAL MODELING 


The jet-induced mixing problem considered in the present study is steady- 
state and incompressible with gravity acting in the vertical negative-x direc- 
tion. The dimensionless forms of the governing equations, with buoyancy force 
resulting from the temperature gradient in the liquid neglected, are 


3u* 


3r*v* 


3x* + r* 3r* = 0 


(7) 


* ij * 


9u*v 

~3? 


3u* 2 

3u‘ 


3p* 

1 


3 2 U* 

3 

fr* ^ 

3x* + 

r* 

3r* = 

3x* 

B Re. 

J 

,3x* 2 

+ r* 3r* 

V 3r* 

3r*v* 2 


3p* 

1 

V* 


1 

3 2 v* 

3 (V 

r* 9r* 

m 

Sr 7 

*' fc i 

2 + 
r*^ 

B 

Re j 

[ 9 x‘ 2 + 

' 3r* V 

3u*T* 


3r*v*T* 


1 

3 2 T* 

3 ( 

r* anV 

3x* 

’ + 

r* 3r* 

B Re^ Pr 

[ 9 x* 2 + 

r* 3r* v 

3rV 


( 8 ) 


(9) 


(10) 


The nondimensional variables in the above equations are defined as 


* r * _ 1 u* _ H_ V* = — „■ 

D’ “O’ u j’ u j’ M 


P - Pq T 

sA T* = 


1 


pu 




Wj 


It Is noted that, for convenience, the gravity term has been subtracted from 
the x-momentum equation by using the static equilibrium equation: 


3p 

_ia 

3x 


-pg 


The relevant parameters in the governing equations are the jet Reynolds number 
(Rej), the Prandtl number (Pr), and tank to jet diameter ratio (B): 


pu^d C p q 

Re j * p ’ Pr " k ' B = d 


(11) 


The following boundary conditions are used to solve the elliptic equations 
(7) to (10). At the centerline, the symmetric conditions are used: 

v * _ o d,y_ _ _ q 

v * °’ 3r* 3r* " U 


Nonslip conditions are applied to the solid walls: 

u* = v* = 0 
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The bottom-wall is assumed to be adiabatic: 



The side-wall is subjected to a uniform heat flux: 



where the wall heat-flux parameter Nh is defined by 


Nh 


q D 
H w 


k<T s - V 


( 12 ) 


The interface Is at saturation temperature T s and is assumed to be wave- 
free and shear-free: 



T* = 1 


at 



The axial velocity at the interface is the condensation-induced velocity and 
is given, from equation (6), by 


u* = 



Pr B Re^ 


at 



(13) 


It should be noted that although buoyancy force is neglected, the velocity and 
temperature fields are still coupled through the condensation-induced velocity 
u c at the Interface. 

For the central jet, uniform velocity and temperature are assumed: 


u* = 1 , v* = 0, T* = 0 

In order to maintain constant liquid level, the volume flow rate of the liquid 
withdrawn from the outer part of the tank bottom should be equal to the jet 
volume flow rate plus the condensation volume flow rate. With the assumption 
of uniform velocity and zero temperature gradient, the boundary conditions at 
the liquid-withdrawn plane are given by 


u* 



1 


+ 



and 


v* = 0, 


ill 

3x* = 


0 
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It is evident that besides Re j , Pr, and B (= D/d), four additional 
parameters associated with the boundary conditions are: 


H 

D* 


Ja, and 
A j 


Nh 


( 14 ) 


The average condensation heat-transfer coefficient, h c , and condensation 
Stanton number, St c , which describe the interfacial heat and mass transports 
are defined as 


IT) h~ 

h. 


c <T s - V 


(15) 


St. = C 


C " pU j C p 


(16) 


By using equation (4) or (5), the average condensation Stanton number can be 
expressed as 


or 


respectively, where 


and 


„0 . 5 


(f£) 


St. = 


Re. B Pr 

J 


(17) 


St. = 


c _ Re. B Pr’ 

J 


(18) 


( Pr . Re J T °ut) . 4 Nh (H 

D\ VD. 


1 + Ja 


0 T out) 


t* _ T out T ,i 

out “ T s - T. 
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From the above analysis, the condensation Stanton number can be expressed 


as 


St, 


f Re 


y 


p r - - 

rr ’ d’ D’ 


out 


Ja, Nh 


(19) 


The effect of Re-j and Pr on St c has been studied in reference 7. This 
paper describes the effects of the other system parameters on the interface 
condensation process. The calculations are conducted for jet Reynolds number 
(Re-s ) equal to 150, 200, and 300 which are In laminar flow range, and for 
Prandtl number (Pr) equal to 1.25 and 2.1 which represent typical values of 
liquid hydrogen and nitrogen, respectively. Also, the parameters D/d, H/D, 
Aout^Aj, Ja, and Nh are varied from 12.5 to 33.3, 0.25 to 1.25, 7.5 to 360, 
0 to 0.2, and 0 to 4, respectively. 


NUMERICAL METHOD OF SOLUTIONS 

The above elliptic partial differential equations are numerically solved 
by a finite-difference method. The finite-difference equations are derived by 
integrating the differential equations over an elementary control volume sur- 
rounding a grid node appropriate for each dependent variable (ref. 8). A 
staggered grid system is used such that the scalar properties, p and T, are 
stored midway between the u and v velocity grid nodes. The bounded skew 
hybrid differencing (BSHD) is Incorporated for the convective terms (ref. 8) 
and the Integrated source terms are linearized. Pressures are obtained from a 
predictor-corrector procedure of the Pressure Implicit Split Operator (PIS0) 
method (ref. 9) which yields the pressure change needed to acquire velocity 
changes to satisfy mass continuity. The governing finite-difference equations 
are solved iteratively by the ADI method with under relaxation until the solu- 
tions are converged. 

Calculations are performed with a nonuniform grid distribution with con- 
centration of the grid nodes in the centerline, near-wall, and near-interface 
regions where the gradients of flow properties are expected to be large. The 
nonuniform grid distribution In axial direction Is generated by using an expo- 
nential function of Roberts' transformation (ref. 10) with the stretching 
parameter equal to 1.02. In the radial direction, the scheme used by Cebecci 
and Smith (ref. 11) with a constant ratio between two adjacent grid spacing is 
used. The 72 by 41 grid nodes which have been shown In reference 7 to give 
reasonable grid-independent solutions are used for all the calculations in the 
present study. Calculations are performed on a CRAY-XMP computer located at 
NASA Lewis Research Center. The convergent solutions are considered to be 
reached when the absolute value of (eqs. (17) - (18))/ ( eq - (17)) Is less than 
0.005 and the maximum of absolute residual sums for each dependent variables is 
less than 10 - ^. 


RESULTS AND ANALYSIS 


Numerical solutions showed that the flowfield near the jet region is 
generally Independent of the parameters Ja, Pr, Nh, and A ou ^/Aj at least 
for their ranges considered in this study. The linear relation between the 
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potential core length to jet diameter ratio (xp/d) and jet Reynolds number 
(Rej) obtained in reference 7, Xp/d = 0.0067 Rej, is confirmed. The present 
calculations further show that tne linear relation holds only when H/d is 
greater than some critical value. As shown in figure 2, larger Rej and 
lower D/d tend to reduce the critical value of H/d. This is because that 
the larger Rej and d/D carry more jet momentum to resist the disturbance 
Induced by the existence of free-surface . 


The axial jet Introduced from the bottom will turn radially outward from 
the centerline as it approaches the interface. Numerical solutions indicate 
that the radial velocity at the interface v* is essentially unaffected by 

Pr, Nh, and A ou ^/Aj. However, the v* Isa strong function of H/D and D/d 

and is slightly affected by Ja. Figure 3 gives the distribution of radial 
Interface velocity (v*) as a function of Ja, H/D, and D/d. The v* is 

increasing with Ja and decreasing with H/D and D/d. The peak value of v* 

is located closer to the centerline when the parameters D/d and Rej are 

larger, and Ja and H/D are smaller. However, it is generally located 
within the region of r* < 0.1. The average value of v* can be described by 
the following equation within +3 percent of the exact numerical solutions: 


v* = 0.2 


Re 


0.5/H' 


j VD 


-0. 25, -0.905 

W 


(1 + 0.226 Ja) 


(20) 


for H/d > 10, Ja < 0.2, 150 < Rej < 600, Nh < 4, and 1.25 < Pr < 2.65. 


The average condensation Stanton number St c as a function of A 0U |-/Aj 
is shown in figure 4 for adiabatic wall condition, Nh - 0. Calculations also 
Include the different outflow locations for the same A ou t/Aj. It_1s 
obvious that under adiabatic (or very low heat leak) conditions, St c is 
essentially Independent of the liquid-withdrawn location and area. _However, 
additional calculations for high heat leak, Nh = 4, Indicate that St c is 
greater if the outflow location is closer to the tank wall. This is because 
the outflow located closer to the wall enhances the liquid circulation of the 
tank and convects more heat from the wall to gain a higher outflow temperature 
at the tank bottom. This information will be very useful in the design of a 
pumping system for the mixer. 

The heat transfer due to condensation is carried away by the radial flow 
motion near the interface. Therefore, the interfacial heat transfer is expec- 
ted to be greater for lower liquid fill level since the radial velocity near 
the Interface wiJJ. be greater. The variation of the average condensation 

Stanton number St c with H/D is shown in figure 5. For an adiabatic wall, 

i.e., Nh = 0, St c is a relatively weak function of H/D. However, with an 
Increase in the wall heat-flux, the average condensation rate decreases sig- 
nificantly as H/D Increases. As for example, figure 5 shows that for Nh = 4 
and Rej = 150, the condensation rate at H/D = 1.25 is reduced by more than 
60 percent compared to that at H/D = 0.25. 

Figure 6 shows the distribution of the Interface temperature gradient for 
different values of H/D and D/d. The interface temperature gradient at the 
central part of the interface, OT*/3x*) s , increases as H/D decreases. 
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However, (3T*/9x*) s decreases with a decrease in H/D at the outer part of 
the Interface. In contrary to the effect of H/D, the interface temperature 
gradient at the central part decreases with a decrease in D/d, but increases 
appreciably at the outer part of the Interface. The average value of the 
interface temperature gradient is greater for smaller D/d because the larger 
jet diameter introduces more jet flow rate to the tank. Also seen from figure 
6, lower H/D and higher D/d result in the condensation being more confined 
in the central part of the interface. Thus, the normal-g (flat interface) 
solutions may be applied to the low-g condition if the liquid fill level is 
low and jet nozzle diameter is small enough. The effect of D/d on St c is 
given in figure 7. For a given H/D, St^/ (d/D)^ approaches an_asymptoti c 
value. Thus, it can be stated that if D/d is large enough, St c will be 
proportional to (d/D) 2 . For Ja << 1 and Nh << 1 , a correlation for predic- 
ting St c (within +5 percent bound) is given by 


St c = 0.977(£ 


-0.025 


Pr 


-0.03 


-2 


( 21 ) 


for 0.25 < H/D < 1.25, 0^85 < Pr < 2.65, 150 < Rej < 600, and D/d > 20. It 
is noted that although St c is Independent of jet Reynolds number Re j for 
Ja << 1 and Nh << 1, this may not be true for larger values of Ja and Nh. 
From equation (21), it is evident that D/d Is a more Important parameter 
affecting the interface condensation than H/D and Pr . If H/D and Pr are 
close to one, St c =■ (d/D) 2 may give a reasonably good prediction of interface 
condensation rate. In that case, St c is nearly a constant determined by the 
system geometry and the Interface condensation rate can be determined by the 
jet volume flow rate. This conlusion also supports the finding in reference 7. 

The_jet subcooling is characterized by Jakob number (Ja). Figure 8 shows 
that St c is linearly increasing with Ja. For an adiabatic wall (Nh = 0), 
the increase of St c with respect to Ja Is very small. The effect of Ja 
on St c is enhanced at larger Nh. However, even for high wall heat flux, 

Nh = 4, the difference in the value of St c between Ja = 0 and 0.2 is still 
less than 4 percent. For most common cryogens such as liquid hydrogen, nitro- 
gen, and oxygen, Ja is generally less than 0.2 for most practical applica- 
tions. Thus, the effect of Ja can be neglected and the work of reference 7 
In which Ja = 0 was_assumed is also justified. It is noted that the explicit 
effect of Ja on St c is directly through the condensation-induced velocity 
u* at the Interface. Thus, the negligible effect of Ja simply means that, 
whether u* at the interface is Included in the calculations or not, the 
effect is negligible. Numerical solutions show that the maximum value of u c , 
which is usually located at or near the centerline. Is generally at least two 
orders of magnitude lower than the jet velocity. For Instance, (u c /uj) max 
Is about 0.0077 for Ja = 0.2, Nh = 0, Rej = 300, Pr = 1.25, H/D = 1.0, and 
D/d = 20. 


Heat transfer through the tank walls tends to generate thermal stratifica 
tion and to Increase the average temperature of the liquid. Thus, an increase 
in the wall heat-flux will result in the reduction of condensation rate. Fig- 
ure 9 shows the linear relation between the average condensation Stanton number 
Si c and the heat leak parameter Nh. The distribution of interface temper- 
ature gradient is shown in figure 10. The Interface temperature gradient 
(8T*/8x*) s is lower everywhere for Nh > 0 (wall heat leak) than that for 
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Nh = 0 (adiabatic wail). It is noted that if Nh is large enough, say 
Nh = 4, a slight vaporization may occur near the tank wall region. 


CONCLUSIONS 

The effects of system parameters on the interface condensation rate in a 
laminar jet-induced mixing tank as shown in figure 1 have been numerically 
investigated. The selected jet velocities were low enough such that the free 
surface was approximately flat. The effect of vapor superheat was assumed to 
be negligible and the interface condensation rate was determined from the 
resulting temperature field in the liquid region. Solutions were obtained by 
solving the nondimensional form of steady-state conservation equations with a 
finite-difference method. Parameters investigated included liquid height to 
tank diameter ratio ranging from 0.25 to 1.25, tank to jet diameter ratio 
ranging from 12.5 to 33.3, liquid outflow to jet area ratio ranging from 7.5 
to 360, and side-wall heat flux parameter varied from 0 to 4, and jet Jakob 
number varied from 0 to 0.2. Under the above stated conditions, the following 
conclusions can be made: 

• The condensation rate is essentially independent of the outflow to jet 
area ratio, A ou t-/Aj, except for high wall heat flux 

• The average condensation Stanton number, St c , linearly Increases with 
Jakob number, Ja, and linearly decreases with the wall heat-flux param- 
eter, Nh. The effect of Ja on St c is enhanced when Nh is greater 
than zero. 

• The average condensation Stanton number, St c , is decreasing with increas- 
ing liquid height to tank diameter ratio, H/D, and tank to jet diameter 
ratio, D/d. The parameters D/d and Nh have stronger effects on St c 
than the parameters H/D and Ja. 

• For Ja << 1 and Nh << 1, St c will be approximately equal to (d/D) 2 if 
D/d > 20 and Pr and H/D are close to one. Thus, the interface conden- 
sation rate can be determined by the jet volume flow rate. 

• Lower values of H/D and d/D will yield condensation more confined at 
the central part of the interface. Thus, normal-g (flat interface) solu- 
tion may be applied to the low-g (curved interface) condition if the 
liquid fill level is low and the jet nozzle diameter is relatively small. 
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Figure 4. - Condensation Stanton number as a function of 
outflow to jet area ratio for Rej = 150, Nh = 0, and Ja = 0 
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Figure 5. - Condensation Stanton number as a function of 
liquid level height for Ja = 0 and Aout^Aj = 357. 
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Figure 6. - Temperature gradient distribution at the interface 
for Rej = 1 50, Pr * 1 .25, Nh = 0, and Ja = 0. 
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Figure 7. - Condensation Stanton number as a function of 
tank to jet diameter ratio for Re j = 1 50, Pr « 1 .25, Nh = 0, 
and Ja = 0. 
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Figure 8. - The effects of Jacob number and wall heat-flux 
parameter on the average condensation Stanton number 
for D/d = 20, H/D = 1 .0, Re j = 300, and Pr = 1 .25. 
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Figure 9. - The effects of wall heat-flux parameter and Jakob 
number on the average condensation Stanton number for 
D/d = 20, H/D = 1.0 , Rej = 300, and Pr - 1.25. 
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Figure 10. - Temperature gradient distribution at the interface 
for H/D = 1 .0, D/d = 20, A out /A i = 357, Re j = 300, and 
Pr * 1 .25. 
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